Localization of a Bose-Einstein condensate vortex in a bichromatic optical lattice 
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By numerical simulation of the time-dependent Gross-Pitaevskii equation we show that a weakly 
interacting or noninteracting Bose-Einstein condensate (BEC) vortex can be localized in a three- 
dimensional bichromatic quasi-periodic optical-lattice (OL) potential generated by the superposition 
of two standing-wave polarized laser beams with incommensurate wavelengths. This is a general- 
ization of the localization of a BEC in a one-dimensional bichromatic OL as studied in a recent 
experiment [Roati et ai, Nature 453, 895 (2008)]. We demonstrate the stability of the localized 
state by considering its time evolution in the form of a stable breathing oscillation in a slightly al- 
tered potential for a large period of time. Finally, we consider the localization of a BEC in a random 
ID potential in the form of several identical repulsive spikes arbitrarily distributed in space. 

PACS numbers: 67.85.Hj,03.75.Lm,03.75.Nt,64.60.Cn 



I. INTRODUCTION 

The possibility of the localization of the electronic wave 
function in a one-dimensional (ID) disordered potential 
as predicted by Anderson in his pioneering work [l[ has 
drawn the attraction of researchers in different areas. 
Different forms of localization have been observed ex- 
perimentally in diverse contexts, such as in electromag- 
netic waves [U, in sound waves H, and also more re- 
cently in quantum matter waves [5|-|8|. Using a cigar- 
shaped noninteracting Bose-Einstein condensate (BEC) 
0, E3] of 87 Rb atoms, Billy et al. Q demonstrated its 
exponential localization when released into a ID waveg- 
uide in the presence of a controlled disorder created by 
a laser speckle. In another experiment, using a cigar- 
shaped noninteracting BEC of 39 K atoms, Roati et al. Q 
demonstrated its localization in a ID bichromatic optical- 
lattice (OL) potential created by the superposition of two 
standing-wave polarized laser beams with different wave- 
lengths. The noninteracting BEC of 39 K atoms was cre- 
ated @ by tuning the inter-atomic scattering length to 
zero near a Feshbach resonance 
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The localization of a BEC in a ID bichromatic quasi- 
periodic OL potential and related to pics have been the 
subject matter of several theoretical [12h25| and exper- 
imental studies. After the pioneering experiments 
on the localization of a ID cigar-shaped BEC, a nat- 
ural extension of this phenomenon would be to achieve 
localization in higher dimensions, e.g., in two (2D) and 
three (3D) dimensions. We address this important issue 
in the present investigation. Using the complete numer- 
ical solution of the Gross-Pitaevskii (GP) equation [26| . 
here we study the localization of a (nonrotating) disk- 
shaped BEC in 2D and also a BEC in 3D, with a small 
nonlinearity, in the presence of bichromatic OL potcn- 
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tials along orthogonal directions. (For zero nonlinearity 
the problem of a stationary state of a c?-dimcnsional — 
d = 2, 3 — nonrotating BEC trivially decouples into d ID 
problems and for a large repulsive nonlinearity the local- 
ization is destroyed [ljll20(.) In 3D another nontrivial 
phenomenon can happen, i.e., the generation of a stable 
vortex state of unit quantum. Here we demonstrate the 
localization of a stable nontrivial vortex state in a 3D 
BEC under the action of bichromatic OL potential along 
the axial direction and harmonic potentials along radial 
directions. We find, as in ID [l^H^l, a repulsive nonlin- 
earity has a strong effect on localization and a not-too- 
large nonlinearity destroys localization in all cases. We 
exhibit results for zero and small nonlinearities. (Effects 
of a weak nonlinearity in Anderson localization have been 
shown experimentally in light waves in photonic crystals 

In the presence of strong disorder, the localized state 
could be quite similar to a localized state of Gaussian 
shape in an infinite potential. However, the more inter- 
esting case of localization is in the presence of a weak 
disorder when the system is localized due to the quasi- 
periodic nature of the potential Q and not due to the 
strength of the lattice. When this happens the localized 
state acquires an exponential tail. The present local- 
ization with an exponential tail in a quasi-periodic OL 
potential with a deterministic weak disorder is a special 
case of Anderson localization in a fully disordered po- 
tential and is well described in the Aubry- Andre model 

In the pioneering study, Anderson considered localiza- 
tion of electron(s) in a ID random potential generated by 
impurity or disorder arbitrarily distributed on a lattice. 
This is of interest in the study of a BEC and we consider 
its localization in a disorder potential in the form of iden- 
tical narrow repulsive spikes (simulating delta functions) 
distributed randomly in space. (This is different from 
a speckle potential where the spikes also have different 
strength and width.) We find that as small as four such 



2 



spikes can localize a noninteracting BEC in ID. This is 
not just of academic interest as such potential can be cre- 
ated in ID and possibly in 2D for a BEC in an atom chip 
[28L |29| . A single repulsive spike separating the BEC into 
two parts is now routinely created [29| . 

There have been theoretical studies on different aspects 
of Anderson localization which are worth mentioning poj . 
Wobst et al. [3(| considered Anderson localization in 
higher dimensions. Sanchez-Palcncia et al. and Clement 
et al. considered Anderson localization in a random po- 
tential [III [It], HH . Damski et al. and Schulte et al. con- 
sidered Anderson localization in disordered OL potential 
[l7| . There have been studies of Anderson localization 
with other types of disorder [32| . Effect of interaction on 
Anderson localization was also studied [HI, [33[ . Ander- 
son localization in BEC under the action of a disordered 
potential in 2D and 3D has also been investigated [HJ|. 

In Sec. [TTl wc present a brief account of the nonlin- 
ear 2D and 3D timc-dcpcndcnt CP equations used in our 
study and of the variational solution of the same under 
appropriate conditions. The generalization of the equa- 
tion to study vortex states in 3D is also presented. In 
Sec. IIIII we present numerical results of localization em- 
ploying time propagation using the semi-implicit Crank- 
Nicolson algorithm. The wave function of the localized 
states has a central Gaussian (variational) form with a 
long exponential tail. First we consider in Sec. MI Al 
the localization of a 2D disk-shaped BEC. We also con- 
sider localization of a 3D BEC and a 3D vortex BEC in 
Dichromatic OL potential(s) in Sec. MI Bl In Sec. MI CI 
we consider the localization of a noninteracting BEC for 
a random potential comprised of arbitrarily distributed 
narrow spikes. In Sec. IIVI we present a brief discussion 
and concluding remarks. 



II. THEORETICAL FORMULATION OF 
LOCALIZATION 

The quasi-periodic bichromatic OL potentials gener- 
ated by two standing-wave polarized laser beams of in- 
commensurate wavelengths in the x direction have the 
following generic forms 



U(x) =^ S i E i COS 2 (fciX), 
i=l 
2 

U(x) = SiEi sin 2 (kix), 



(1) 

(2) 



where Sj, i = 1,2, are the amplitudes of the OL potentials 
in units of respective recoil energies Ei = 2n 2 Ti 1 / (m\ 2 ), 
and hi = 27r/ Aj, i = 1, 2 are the respective wave numbers 
and Aj are the wave lengths, h(= h/2ir) is the reduced 
Planck constant, and m the mass of an atom. 

If we have a single periodic potential of forms (TT|) and 
@ with S2 = 0, the solution of the Schrbdinger equa- 
tion cannot be localized. One can have localization if a 



second periodic component with a different frequency is 
introduced in Eqs. (fTJ) and @. These localized states 
are not the gap solitons, which are localized states in 
the solution of a nonlinear Schrbdinger equation with a 
repulsive nonlinearity appearing in the band-gap of the 
spectrum of the linear Schrbdinger equation [35| . 
The BEC in 3D is described by the CP equation 
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where g = 4nh 2 aN/m, j |</>(r, r)\ 2 dr = 1, r the time, N 
the number of atoms, V(r) is the trap, and a is the atomic 
scattering length. With three quasi-periodic, bichromatic 
OL potential in x, y, and z directions, after canceling the 
factor h 2 /m from both sides of Eq. ([3]), the GP equation 
in explicit notation becomes 



■ d<j)(x,y,z,t) 
dt 



(d 2 + d 2 +d 2 )+V(x) + V(y) 



V{z)+g\4>{x,y,z,t)f 



where g = AiraN, d x 's denote space derivatives, and 
time t = rh/m is now expressed in units of m/h. 
Note that Eq. ([!]) is not expressed in dimensionlcss 
units. The variables x, y, z, Aj are in actual units of 
length (X), \4>\ 2 is in units of L~ 3 with normalization 
I-00 I-00 I 00 \<K X > V> z > t)?dxdydz = 1. InEq. ©the 
scaled potentials V(x) = U(x)m/h 2 are now defined by 
one of the two following expressions: 
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In case of a 2D and 3D BEC, instead of having the same 
potential, e.g., (J3J) or (J5J) in different directions, one can 
choose different potentials along different directions. 

For axially-symmetric traps, the GP equation can be 
easily generalized to include a vortex state, as shown in 
Ref. floL l36l| . in axially-symmetric coordinates r = (p, z), 
where p is the radial coordinate and z the axial coordi- 
nate. To obtain a singly quantized vortex state of angular 
momentum Ti around z axis, one has to explicitly intro- 
duce a phase (equal to the azimuthal angle) in the wave 
function. (Vortex states of higher angular momentum are 
unstable and decays into multiple states of angular mo- 
mentum Ti.) This procedure introduces a centrifu gal term 
in the GP equation representing a vortex state |l0l l3(| . 
Thus we can study a localized quantized BEC vortex of 
unit angular momentum in the axially-symmetric poten- 
tial, where the bichromatic OL potential V(z) is placed 
along the axial z direction and an harmonic trap p 2 /2 
along the transverse radial p direction. The modified GP 
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equation for such a vortex is given by [3] 
,d(f>(p,z,t) 
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where we have explicitly included the angular momentum 
dependent centrifugal term l/(2p 2 ) [I?]]- Because of the 
centrifugal term, the density of the vortex should be zero 
along the z axis. 

If there is a strong harmonic trap in the axial z direc- 
tion, one can derive a reduced quasi-2D GP equation for 
a disk-shaped BEC, which can be written in dimension- 
less harmonic oscillator units as 1381 



■ d4>{x,y,t) 
1 dt 



g\(j)(x,y,t)\' 



d 2 



, + V(x) + V(y) 
<t>fay,t), (8) 



where g = 2\ / 2iraN and V(x) is given by Eqs. 
(J5J or ([BJ . In Eq. © the normalization is 
J Too JT 10(^1 Vi t)\ 2 dxdy = 1. In Eq. lengths are ex- 
pressed in units of l z = y/h/moj z with lu z the frequency 
in z direction, (j){x, y, t) in units of l~ 2 and time in units 
of u)- 1 . 

Finally, for the sake of completeness we note that if 
there is a strong harmonic trap in transverse y and z di- 
rections, one can derive a reduced quasi-lD GP equation 
for a cigar-shaped BEC, which can be written in dimen- 
sionlcss harmonic oscillator units as 1381 



.d(j)(x, t) 
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where g — 2aN and V(x) is given by Eqs. (JS|) or ©. 
Now the normalization is J_ ca \(j>(x,t)\ 2 dx = 1. Equation 
© has been used [TH, [2(| for the study of localization in 
ID in bichromatic OL potential. In Eq. © lengths are 
expressed in units of I = ^Jh/muj±_ with io± the frequency 
in transverse y and z directions, <fi(x,t) in units of I -1 and 
time in units of . 

Although we use potentials (J5J) and (JB]) in our study, 
there is some difference between these two potentials. Po- 
tential © generates a different type of localized states 
compared to potential ([5]). Potential ([6J has a local min- 
imum at the center, consequently stationary solutions 
with this potential have a maximum there. However, 
potential ([SJ has a local maximum at the center corre- 
sponding to a minimum of the stationary solution. 

Usually the stationary localized states formed with 
quasi-periodic OL potentials (J5j and (JHJ) occupy many 
sites of the quasi-periodic OL potential and have many 
local maxima and minima. For certain values of the pa- 
rameters, potential (jBJ leads to localized states confined 
practically to the central cell of the quasi-periodic OL po- 
tential. When this happens, a variational approximation 



with Gaussian ansatz leads to a reasonable prediction for 
the localized state in the central region. (However, it has 
an exponential tail at large distances.) 

To derive a simple variational solution of linear Eqs. 
(g]), ©, and (0 in 3D, 2D, and ID, respectively, in a 
unified fashion, we adopt the convenient notation r = 
(x,y,z) in 3D, = (x,y) in 2D, and = (x) in ID. The 
stationary form of the linear Schrodingcr equations Q , 
(HJ, and (|9]) (with id/dt replaced by a chemical potential 
p) with potential ([BJ can be derived from the following 
Lagrangian 



L = 



/#(r)| 2 --|V0(r)| 2 



V>Mr)| 2 - f |<Xr)r 



dr — p, 



(10) 



by demanding SL/54> — SL/5p — 0. To apply the varia- 
tional approximation we use the Gaussian ansatz [39| 
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where d is the dimension of space, x\ = x, X2 = y, x% = z, 
and the variational parameters are the norm j\f, width w, 
and p. This ansatz implies that the center of the station- 
ary state is placed at the local minimum at Xi = 0, i = 
1, 2, 3 of the quasi-periodic OL potential. The substitu- 
tion of ansatz (fTTj) in Lagrangian (fTU)) leads to 

L = p{N - 1) - + dj\f^2 ^-[exp(-a 2 U! 2 ) - 1] 
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where Ai = 2n 2 s i /\ 2 ,a i = 27r/Aj. The first variational 
equation from Eq. (fT2"]) . dL/dp = 0, yields j\f = 1, which 
will be used in other variational equations. The second 
variational equation dL/dw = 0, yields 



1 = 2u 2 A iW 4 exp(-a 2 w 2 ) 
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and determines the width w. The last variational equa- 
tion dL/dN = 0, yields 



^-£^[« P (-aM-l] + 9 y 

l—l 
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which determines the chemical potential. 

In addition to considering a disorder potential in the 
form of a bichromatic lattice, we also consider the follow- 
ing disorder potential in the form of randomly distributed 
S repulsive identical Gaussian spikes in ID along the x 



V(x) =^B 8 exp[-c i (x~/3 J ) 2 



(15) 
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where Bi is the amplitude of the Gaussian spike, c, is its 
width, and /?, is its random position. If the spikes are 
placed at constant periodic spacing, no localization can 
be obtained. There will be localization if the positions 
Pi are random. In 2D and 3D the appropriate potentials 
are V(x) + V(y) and V(x) + V(y) + V(z). We study the 
localization of a BEC in random potential (|15|) . 

For an analytical understanding of the problem, next 
we present a variational analysis. As the potential (|T5|) is 
Gaussian a variational analysis based on the Gausssian 
ansatz is fully intcgrablc. Only the potential term in the 
Lagrangian gets modified and in this case the third term 
in the Lagrangian of Eq. (|T2"j) becomes 
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FIG. 1: (Color online) Numerical density 4> 2 (x) (line) and its 
exponential fit (symbol) vs. x for Ai = 15, 10 and 5. 



and Eq. (|T3| for width gets modified into 
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where ji = 1 + qui 2 . The expressiooafor chemical poten- 
tial becomes 
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III. NUMERICAL RESULTS 

We performed numerical simulation employing 
imaginary- and real-time pr opag ations with Crank- 
Nicolson discretization scheme [U El using adequately 
small space and time steps necessary for obtaining 
converged solutions. In practice we used space and time 
steps smaller than 0.025 and 0.0005, respectively, and 
sometimes as small as 0.0025 and 0.00002, respectively. 
We use the FORTRAN programs provided in Ref. [38[ 
for our purpose. For checking the consistency of our 
calculation we compared our real-time results with 
imaginary-time results and we verified that the two sets 
of results were in agreement with each other. Because of 
the oscillating nature of the bichromatic OL potential, 
great care was needed to obtain a localized state 
precisely. The accuracy of the numerical simulation was 
tested by varying the space and time steps as well as 
the total number of space and time steps. Although we 
used the time-dependent GP equation for the study of 
localization, all results reported in this paper, except 
those in Fig. |H are stationary results independent of 
time. In Fig. U we study time-dependent stability 
dynamics of the localized states obtained with real-time 
propagation. 
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FIG. 2: (Color online) (a) Numerical and (b) variational den- 
sity 4> 2 (x, y) vs. x and y from a solution of Eqs. (JSJ and (0 for 
a disk-shaped BEC for g — 0. (c) Numerical density 4> 2 (x,y) 
vs. x and y and (d) its contour plot from a solution of Eqs. 
(IH) and ((6]) for a disk-shaped BEC for g = 2. (e) Numerical 
density (f> 2 (x,y) vs. x and y and (f) its contour plot from a 
solution of Eqs. (JHJ and © for a disk-shaped BEC for g — 2. 
Quantities (j) 2 (x,y) and x,y are all in dimensionless harmonic 
oscillator units. 
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A. 2D Bichromatic Optical Lattice 

To study the localization of a BEC in 2D and 3D with 
potentials ([5]) and ([5]), we set the ratio A 2 /Ai = 0.86 
(roughly the same ratio A2/A1 as in the experiment of 
Roati et al. Q). First we consider the solution of 
Eq. ((HJ for a disk-shaped BEC. To understand the na- 
ture of these localized states, we consider the localized 
states with larger values of Ai. Such states with a large 
S2I s\{= 1) occupy a small number of OL sites and hence 
their numerical simulation can be performed relatively 
easily All results reported in Sees. IIII Al and lHI Bl are ob- 
tained with the following parameters in the bichromatic 
OL potentials §5§ and ©: Ai = 5, A2/A1 = 0.86, and 
Si = s 2 = 4. It now remains to be seen if with this 
set of parameters the localized state is in the limit of 
weak disorder with an exponential tail. In ID this would 
mean <fi(x) ~ exp(— x/L\ oc ), where L\ oc is the localization 
length. In Fig. Q] we plotted the numerical probability 
<fi 2 (x) vs. x on a log scale, together with a fitting ex- 
ponential function with L\ oc = 0.75, 1.5 and 2.2, respec- 
tively, for Ai = 5, 10 and 15. These localization lengths 
are large compared to the root mean square (rms) size 
of the localized states, which are 0.53, 1.05, and 1.55 for 
Ai = 5, 10 and 15. The wave functions shown in Fi g. [1] 
are identical to those shown in Fig. 1 (a) of Rcf. [20[ 
where the central part of these wave functions are well 
fitted to Gaussian variational approximations. 

In Figs. [2] (a) and (b) we plot the results for den- 
sity (j) 2 (x,y) vs. x and y from numerical and variational 
solutions of Eqs. (|8) and ^ for g = 0, respectively. 
(The g = case is trivial as the 2D and 3D Eqs. (|SJ) 
and (g]) decouple into two and three ID equations, re- 
spectively. Nevertheless, this case allows us to test care- 
fully the numerical programs by comparing the varia- 
tional result with the two numerical results obtained by 
real- and imaginary-time propagation.) The variational 
width (w = 0.4864) obtained from a solution of Eq. ([13")) 
produced the density in good agreement with the numer- 
ical density in Fig. [5] (a). In Fig. H](c) we plot the results 
for density (fi 2 (x,y) vs. x and y from a numerical solu- 
tion of Eqs. (JSJ and (JB|) for the same set of parameters 
as in Fig. (a) but with g = 2. (As in ID [H Soj], re- 
pulsive nonlinearity g has a strong effect on localization. 
As g is increased to a small positive value, the localized 
states occupy more and more lattice sites with larger spa- 
tial extension.) In Fig. [2] (d) we show a contour plot of 
the density shown in Fig. [5] (c) . Because of the nonlin- 
ear inter-atomic repulsion, the essentially single peak of 
Fig. [U (a) now transforms into a multi-peak structure, 
although the central peak is still the strongest. Finally, 
in Figs. [5](e) and (f) we plot </> 2 (cc, y) vs. x and y and its 
contour plot from a numerical solution of Eqs. © and 
([5]) for g = 2. Due to a change in the nature of the poten- 
tial from Eq. ([6]) (sine, minimum of potential at origin) 
to Eq. ([S]) (cosine, maximum of potential at origin) the 
central region has a minimum of density in this case and 
not a maximum as in Figs. [5] (a) and (c). The difference 
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FIG. 3: (Color online) (a) Density (j) 2 (x,y) of a disk-shaped 
BEC from a numerical solution of Eq. J5} with (b) its contour 
plot for potential ([6]) along x direction and potential ([5]) along 
y direction for g — 0. (c) Same as (a) but with g = 2. (d) 
Same as (b) but with g — 2. All quantities are in oscillator 
units. 
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FIG. 4: (Color online) Time evolution of the rms sizes of 
the BECs depicted in Figs. [2] (c) and (e) as the potential is 
multiplied suddenly by a factor of 1.02 or 1.05. All quantities 
are in oscillator units. 



in the structure of the BEC densities can clearly be seen 
from the respective contour plots. 

We also calculated the chemical potential of the states 
illustrated in Figs. [2] The numerical result for energy 
for potential ^ of Fig. [2] (a) is 4.747, to be compared 
with the variational result of Eq. ([T4l 4.791, calculated 
with width w = 0.4864 obtained by solving Eq. (|13l) . 
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This agreement between the numerical and variational 
results of the BEC density for potential © in Fig. H](a), 
and of the respective energies, provides assurance about 
the accuracy of the numerical code used in simulation 
in our investigation. We also calculated the (numerical) 
chemical potential of the BEC density displayed in Figs. 
|2](c) (fi = 5.102) and (e) (p = 5.103). The two chemical 
potentials arc practically the same. A similar finding was 
noted in the ID case as the potential was changed from 
sine to cosine type [2(| ■ 

As we are considering two potentials ([5]) and ([6]) in 2D 
and 3D it is possible to have a distinct potential along 
each axis. We consider this possibility in the case of 
the localized state in a disk-shaped BEC where we take 
potential ([5]) along x direction and potential along y 
direction with g = 0. The density and the contour plot 
of the BEC in this case are shown in Figs. [3] (a) and (b). 
Just to illustrate the effect of a small nonlinearity on the 
localized state we repeated the numerical simulation in 
this case with a nonlinearity g = 2. The density and 
the corresponding contour plot of the resulting localized 
state are shown in Figs. [3] (c) and (d). The density in 
Fig. [3] (c) is quite similar to that in Fig. [3] (a), with 
the exception that the density in Fig. [3] (c) extends over 
more sites due to the repulsive nonlinearity. 

Next we explicitly demonstrate that the densities 
shown in Figs. [2] corresponds to stable states. For this 
we consider the time evolution of the BECs illustrated in 
Figs. [2](c) and (e) in a slightly altered potential, e.g., the 
potential obtained by multiplying the original potentials 
by a factor of 1.02 and 1.05. The resultant rms sizes as 
calculated by the real-time propagation routine are plot- 
ted in Fig. [U The curves in Fig. [4] are labeled by Fig. [2] 
(c) or (e) and the factor 1.02 or 1.05 which multiplies the 
potential. In all cases the rms sizes execute breathing os- 
cillation over a long period of time as shown in Fig. [4]and 
this demonstrates the stability of the localized states. 



B. 3D Bichromatic Optical Lattice 

Now we consider a few cases of the localized states 
in 3D as obtained from a solution of Eq. ^ together 
with potential ([5]) or © along three orthogonal direc- 
tions. First, as in 2D, we consider the solution of Eq. 
(0| with potential © for g = 0. The result is illus- 
trated in Figs. [5] (a) and (b). In Fig. [5] (a) we plot 
three sections of the density — $ 2 (.x) = </> 2 (x,0, 0) vs. 
x, $ 2 (r) ee (j> 2 (r/V3,r/V3,r/V3) vs. r, and $ 2 (p) ee 
</> 2 (0, p/y/2, p/y/2) vs. p — together with the variational 
result. Here $ 2 (x) corresponds to the density in the axial 
x direction with polar angle 9 = 0, <& 2 {r) that in the di- 
agonal direction with polar angle 6 = 7r/4 and azimuthal 
angle <p = it/ 4, and $ 2 (/o) that in the transverse direction 
with polar angle 9 = 7r/2 and azimuthal angle ip = 7r/4. 
In Fig. [5] (b) we show the 3D contour plot (obtained us- 
ing Mathcmatica) of the BEC density showing the actual 
shape of the localized state. (The value of the density at 



1.6 




FIG. 5: (Color online) (a) Sections of density <j> 2 (x, V, z ) 
for a 3D BEC from a solution of Eq. Q with the (b) 
3D contour plot of density <p 2 (x,y,z) for potential ([6]) with 
g = 0; (c) Same as (a) with g — 2; (d) Same as (b) with 
g — 2; (e) Same as (a) with potential ([5]) with the parameters 
Ai = 5,A 2 /Ai = 0.86, si = s 2 = 4 and g = 2; (f) Con- 
tour plot for (e). The sections plotted in the left panel are 
® 2 (x) = 2 (x,O,O) vs. x, $ 2 (/5) = cf> 2 (0,p/V2,p/V2) vs. p, 
and $ 2 (r) = 2 (r/V3, r/y/3, r/i/3) vs. r. For potential © 
the variational result for <J? 2 (r) vs. r is also plotted. 



the boundary of the plot is 0.001 in Figs. [S] (b), (d) and 
(f).) As g = in Fig. [5] (a), the 3D wave function is 
trivial and decouples in the three directions. The varia- 
tional result in this case is in good agreement with the 
numerical result for the density in the diagonal direction 
— $ 2 (r). Next we consider the nontrivial 3D case with 
g = 2 for potential In this case the 3D solution has 
no ID counterpart. In Figs. [S] (c) and (d) we plot the 
sections of the densities and the 3D contour plot, respec- 
tively, for potential ([6]) with g = 2. The shape of the 
BEC density is similar to the g = case with a maxi- 
mum at the origin, the only difference being that in Fig. 
[S] (c) the density extends over a larger distance in space 
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FIG. 6: (Color online) (a) 2D contour plot of the BEC density 
(f> 2 (p, z) for a localized 3D vortex in Dichromatic OL potential 
([6]) with g = 5, obtained by solving Eq. Q. The vortex core 
appears as the line p — 0. (b) Same as (a) in bichromatic OL 
potential ([5|. As potential ([5]) has a maximum at z = the 
matter density is zero there. 



due to the repulsion introduced by a positive g value. 
The density has secondary maxima in adjacent OL sites 
in this case. Finally, in Fig. [S] (e) and (f) we show the 
results for the density and its 3D contour plot, respec- 
tively, with potential §5§ and g = 2. Now the density 
has a minimum at the origin in contrast to the maxima 
in Figs. [5] (a) and (c). Also, in Figs. [5] (e) and (f) the 
density is zero along the three axes (explicitly shown in 
x direction with <J> 2 (x) = 0). 

Next we consider the localization in a 3D bichromatic 
OL potential of a noninteracting BEC vortex with g = 5 
rotating around the z direction with unit angular mo- 
mentum from a numerical solution of Eq. (|7|). In this 
case we consider potentials (JB|) and along the z direc- 
tion. For a localized 3D vortex state rotating around the 
z direction and described by Eq. (0 under the action of 
the bichromatic OL potentials ([6j and ((5j) , a vortex core 
of zero density passing through the origin develops along 
the z axis. To illustrate the vortex state we solve Eq. ([7]) 
with potentials © and ([5]) for g = 5. In Figs. [B] (a) and 
(b) we show the contour plot of the density 4>(p, z) for 
potentials © and ([5]), respectively. For potential (O, m 
addition to the vortex core along the z axis, the density 
is also zero along the z = line. 



C. Random Spike Potential 

Now we consider the localization of a BEC in the ID 
random potential (|15[) . For small values of the ampli- 
tude Bi, a large number S of spikes is needed for a good 
localization. However we consider the minimum number 
S(= 4) of spikes for localization. For a robust localization 
with S = 4, we take Bi = 8 and Cj = 100 corresponding 
to a Gaussian spike of very small width. Next we had 
to choose the random positions ft. Different sets of un- 
evenly distributed ft produced localization and here we 
take ft = -12,-6,6,14 for i = 1,2,3,4. This choice of 
the points sets the center of the localized BEC approx- 
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FIG. 7: (Color online) The numerical (num) and variational 
(var) BEC density 4> 2 ( x ) vs - x f° r the random potential (|15|) 
with Bi — S,Ci — 100, Xio = —12, X20 = —6, a-30 = 6, £40 = 14 
and g — 0. The potential V(x) is also plotted in arbitrary 
units. 

imately at x = 0. In Fig. [7]wc plot the localized BEC 
density 4> 2 (x) vs. x for nonlinearity g = 0. The varia- 
tional and numerical chemical potentials are 0.03739 and 
0.03093, respectively. The plot of the potential (in arbi- 
trary units) is also shown in Fig. [7J With this potential 
the localization is destroyed for a very small nonlinearity 
g{p~> 1), unless the number of spikes is increased and we 
do not consider the localization of an interacting BEC 
here. 



IV. SUMMARY 

In this paper, using the numerical solution [38j of the 
GP equation [26| in 2D and 3D, we studied the local- 
ization [H, Q of a noninteracting and weakly interacting 
BEC in a quasi-periodic bichromatic OL potential along 
different axes. We considered two analytical forms (sine 
and cosine) of the OL potential, e.g., (O and © and 
considered the same or different potentials along different 
axes. First we consider the localization of a disk-shaped 
BEC with a small repulsive nonlinearity under the action 
of the bichromatic OL potential J5j or J6J. We consid- 
ered the localization of a full 3D BEC under the action 
of bichromatic OL potentials (|5|) or $Q along the three 
axes with or without a small nonlinear atomic interac- 
tion. The increase of nonlinearity destroys localization 
in all cases considered fl9l . |2(I l24j . We also considered 
the localization of a noninteracting BEC vortex of unit 
angular momentum with a 3D bichromatic OL potential 
(J6J) or (J5j) along the axial direction and harmonic traps 
along the transverse radial directions. The clear stable 
vortex cores in 3D as shown in Figs. [5] (a) and (b) are 
the most important findings of this paper. Finally, we 
showed that one can have localized noninteracting BEC 
states under the action of random potentials taken in the 
form of repulsive spikes randomly distributed in space. 
All these localized states are found to be dynamically 
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stable. 

Wc hope that the present work will motivate new stud- 
ies, specially experimental ones, on localization in a more 
realistic 2D and 3D BEC under the action of bichromatic 
OL potentials along different axes. Specially challenging 
is the localization of a singly-quantized vortex BEC in 
3D as predicted in this work. The consideration of a 
vortex lattice of BEC [4l[ in bichromatic OL potentials 
is of great interest also. It remains to be seen if one 
has vortex pinning in a rotating BEC as observed in a 
monochromatic OL potential superposed on a harmonic 
trap [42|. The effect of a bichromatic OL potential on a 
two-species mixture of BECs [43| is also of interest. In 
the bichromatic OL potential considered in this paper, 
the disorder term is taken in the form of a periodic po- 
tential with an incommensurate wave length. We also 



considered a disordered potential in the form of narrow 
identical repulsive spikes distributed randomly in space. 
Further work need be done to study the effect of differ- 
ent types disorder terms. The 2D and 3D localization of 
BEC under the action of bichromatic OL potentials and 
of a random potential — as predicted in this paper — 
may find application in other contexts in 2D and 3D, for 
example, in the localization of electron waves in 2D and 
3D lattice as well as in the localization of electromagnetic 
waves, or sound waves in the presence of disorder. 
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